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We consider the gluon propagator D(p 2 ) at various lattice sizes and spacings in the case of pure 
SU(3) Yang-Mills gauge theories using the Landau gauge fixing. We discuss a class of fits in the 
infrared region in order to (in)validate the tree level analytical prediction in terms of the (Refined) 
Gribov-Zwanziger framework. It turns out that an important role is played by the presence of 
the widely studied dimension two gluon condensate (A 2 ). Including this effect allows to obtain an 
acceptable fit up to 1 a 1.5 GeV, while corroborating the Refined Gribov-Zwanziger prediction for 
the gluon propagator. We also discuss the infinite volume extrapolation, leading to the estimate 
D(0) = 8.3 ± 0.5 GeV -2 . As a byproduct, we can also provide the prediction (g 2 A 2 ) « 3 GeV 2 
, obtained at the renormalization scale a = 10 GeV. 
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Although confinement of gluons in pure Yang-Mills gauge theories should be understood in a gauge invariant 
setting, one should also observe in some way the messengers of confinement in gauge variant quantities. In particular, 
let us assume that we have fixed our gauge freedom by means of the Landau gauge, d^A a ^ — 0. We expect that the 
nonpcrturbative physics responsible for confinement will also reveal their influence on the n-point correlation functions 
of the gauge fixed theory. The most elementary, albeit already nontrivial, example of such correlation function is the 
l— ~~ '■ gluon propagator. Due to the transverse nature of the Landau gauge, we can write 

^> ; (A«(p)Al(-p)) = D{p 2 )5 ab (V - , (1) 

| and focus attention on the scalar quantity D(p 2 ). This quantity has been the topic of a lot of investigations in the 
. recent past, be it from numerical [H4lH or analytical viewpoint [l6r{37l |. All data and analytical estimates agree on 
the fact that the gluon is infrared suppressed. There is still some discussion whether it actually vanishes at zero 
momentum or not, but most lattice data seems to indicate that it does not. 
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An important asset in the computation D{p 2 ) is the role played by Gribov copies jig, [38J. In principle, the gauge 
configuration is in the (absolute) Landau gauge if it corresponds to the absolute minimum of the functional 



R[A] = A 2 min = mm / d d x{Alf . (2) 

J_l ■ uGSU(N)/ 



The set of absolute minima defines A, the so-called Fundamental Modular Region (FMR). It is then an easy exercise 
to show that A^ is part of the Gribov region, defined as 

n = { A; d^Al = 0, M ab > } , (3) 

with Ai ab the (Hermitian) Faddeev-Popov operator, defined by 

M ab = -d a D ab = -(S ab d 2 + g f abc A c ^) . (4) 

Notice that the requirement AA ab > is necessary to remove many redundant gauge configurations, as the 
transversality condition d^A^ = has multiple solutions along each gauge orbit. Said otherwise, the Landau gauge is 
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plagued by the existence of Gribov copies. It is important to mention that S ft does not necessarily mean that 
corresponds to the absolute minimum of i?L4]; it can also constitute a relative minimum. Said otherwise, the FMR A 
is a subset of the Gribov region ft. This means that the Gribov region f2 still contains gauge copies, see for instance 
[39M4ll | . It can be shown that Q is convex, bounded in all directions and that is crossed by any gauge orbit [2(J EH, IHJ ■ 



In lattice computations, one fixes the Landau gauge numerically by searching for the "best" solution of the 
minimization procedure for R[A\. As such, one hopes to bring each configuration as close as possible to a gauge 
equivalent one in A. 

In the continuum, it appears to be an incredibly difficult task to implement the absolute Landau gauge. In a first 
approximation, one uses the perturbative Faddeev-Popov action, 



^YM+gf 



\ J d d xF^ + J d d x (b a d^+^d fl Dfc b ) , (5) 

which just implements d^An — 0, by means of the equation of motion of the auxiliary 6-field. The corresponding 
Jacobian determinant is represented by the ghost term c a d fl D'^ b c b . This approach completely ignores the existence 
of Gribov copies, but it is perfectly well-suited for the perturbative quantization of gauge theories. 

A second approximation involves the restriction of the allowed gauge configurations to the Gribov region f2, which 
already partially resolves the issue of gauge copies. It was worked out at lowest order in a saddle point approximation 
in [ll| and later on generalized to all orders in pjKUj]. We shall skip the details, and just mention the result, i.e. the 
local action 

SgZ — SyM+gf (6) 
i / jd^f—acci { T^&b , „bc\ --flCQ / r^ab , bc\ ^ rabco — a/e r^bd d, „ck 2 „ rabc Aa ( , J)c , — bc\ j ( j\r2 -i \ „,4\ 

+ / d x [ip v df, [Dp <p v )-U3 v (D^ u v ) - gj d^ui v D ^ c ip v -7 gj A^ + tp^ J - d [M - 1) 7 j , 
which contains the (Gribov) mass parameter j 2 . This parameter is not free, but self-consistently fixed by means 
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of the so-called horizon condition [lj], El], which reads in its local version (gf abc A^(ip b ^ + Tp b ^)) = -2d(N 2 - l)-y 2 . 
Upon solving, the horizon condition shall give j 2 oc Aqcd- A crucial feature of this local formulation of the 
restriction is that we can control its ultraviolet behaviour, i.e. the action defines a renormalizable quantum field 
theory 0, HI, IH, S3 • As such, a consistent computational framework is obtained. So far, no one has been able to 
improve upon this restriction, in particular it is unclear whether it would be possible to implement the restriction to 
the FMR A to completely overcome the gauge fixing ambiguity. We can only refer to the conjecture of stating 
that at the level of expectation values, no difference will be found upon restricting to £1 or to its subset A. Anyhow, 
the presence of the mass parameter j 2 oc Aqcd will clearly be generating nonperturbative effects in e. g. gluon and 
ghost propagator. 

A well-known important source of nonperturbative effects in gauge theories are condensates, viz. vacuum ex- 
pectation values of certain local operators. Next to the famous gauge invariant condensate (F 2 U ), of paramount 
importance for phenomenological applications (45j . recent years 1 have also witnessed an increased interest in the 
dimension two condensate (A 2 ) in the Landau gauge [H, [50], and related to it the issue of 1/Q 2 power corrections 
48]. The latter corrections would correspond to an extension of the usual SVZ sum rule study of physical 
correlators. Some important early contributions to this field of research can be found in, for example, [5ll459| . 
These works were based on renormalon analyses, lattice considerations of the interquark potential and condensates, 
nonperturbative short distance physics, .... Also at the propagator level such power corrections were identified in [60| . 

From the definition ([2]), it is clear that (^4^^) is a gauge invariant quantity by construction. This leads very 
naturally to the introduction of (A 2 ) in the Landau gauge since we can write [6 If 



A 2 - - I /,"' .* 

^min 2 ' 



0(A 4 ), (7) 



from which it easily follows tha t (A 2 ^- ) = (A 2 ) in the Landau gauge. This condensate then made its appearance in 
a variety of works, see e.g. [22I [63 - l79j | . In the works [4£|[5(|, the relation was explored between this condensate and 



1 The d = 2 gluon condensate was already considered in [46l , [47| . 
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magnetic degrees of freedom, which are generally believed to play an important role for confinement. Recently, this 
was further investigated by looking at the electric and magnetic components of (A 2 ) at finite temperature, hinting 
towards an interesting connection with the phase diagram [70]. 

Measurements of ( A 2 ) at T = have been obtained using the lattice gluon propagator and the Operator Product 
Expansion (OPE) in [gj], based on earlier work [62|, giving the following estimate 

(g 2 A 2 ) = 5.ll?;J GeV 2 (8) 

at the renormalization scale fi — 10 GeV. {A 2 ) also appeared as a source of power corrections in e. g. [ 78l. [79| . An 
independent estimate using the OPE and the quark propagator in a quenched lattice simulation gave [77j 

(g 2 A 2 ) = 4.4 ± 0.4 GeV 2 . (9) 

An ab initio calculation of (A 2 ) was presented in [gH It was shown that it is possible to construct an effective 
potential for (A 2 ) which is consistent with the renormalization (group) (6f| [66|. A nonvanishing condensate due 
to dimensional transmutation was favoured as it lowered the vacuum energy. Using a resummation of Feynman 



diagrams, more evidence for (A 2 ) ^ was given in 67 1 



The extension of the effective potential formalism to the Gribov-Zwanziger case was first tackled in [22|. More 
recently, it also became clear that other d = 2 condensates can play an important role in the Gribov-Zwanziger 
formalism. When the dynamics of the extra fields is taken into account, next to (A 2 ) other dimension two condensates 
appear quite naturally (23l . [24| , and these condensates alter the behaviour of the propagators quite drastically. In 
particular, the d = 2 condensates related to the auxiliary fields {<^ b ,^ b , wj 6 } give a ghost propagator behaving 
like ~ 1/p 2 in the infrared, while the gluon propagator is suppressed and tends to a nonzero constant at very low 
momentum. This framework is now known as the Refined Gribov-Zwanziger (RGZ) formalism, which is a dynamical 
improvement of the original Gribov-Zwanziger approach. In [23|, HH], the effects of the condensate (<^V£f — uJ^ b w^ h ) 
were explored by means of variational perturbation theory. (A 2 ) was left out of this analysis for simplicity, as the 
qualitative conclusions about the deep infrared behaviour were not depending on this condensate, but it was already 
discussed in [24[ that in principle it can be included. In [go], a more complete treatment will be presented. We shall 
not dwell upon details here, but focus on the form of the tree level propagator in the presence of these condensates, 
which is 

2 p 2 + M 2 

(P ' = p 4 + (M 2 + m 2 ) p 2 + 2g 2 N 1 i + M 2 m 2 ' (1 ' 

where M 2 is the mass scale related to the d — 2 condensates in { </?° b , Tp°^ , w" 6 , ZtJ° } , m 2 to (A 2 ) and 7 4 is the Gribov 
parameter. We shall introduce the shorthand A 4 = 2g 2 Nj 4 + M 2 m 2 . 



The aim of this paper is to find out whether the propagator (|10j) can reproduce not only qualitatively the gluon 
propagator, but that it also works out well at the quantitative level. We shall therefore analyze the lattice gluon 
propagator in pure SU(3) Yang-Mills gauge theories and investigate to what extent the propagator (fTU)) can match 
the data, by treating the mass scales m 2 , M 2 and 7 4 as fitting parameters. 

The paper is organized as follows. In Section II, we summarize the introduction of the lattice gluon propagator and 
discuss its renormalization. In Section III, we analyze the propagator and discuss a class of fits related to the RGZ 
propagator (fTQ| . It shall turn out that none of the parameters m 2 , M 2 or 7 4 can be put equal to zero to find a decent 
fit, which thereby shows that the analytical (Refined) Gribov-Zwanziger restriction is in compliance with the lattice 
data, but that it needs to be complemented with the condensate (A 2 ), as well as the RGZ mass scale M 2 , related to a 
new d — 2 condensate in the {<^ b ,^f , fields |24| , which is crucial to find D(0) > 0. We shall also derive an 

estimated value for the condensate (g 2 A 2 ), and we shall see that it compares acceptably well with the OPE estimates 
(jl]) and (jSJ). In addition, we can also derive an infrared gluon mass scale, which lies in the same ballpark as other 
values in the literature. We end with conclusions in Section IV. 



II. THE LATTICE GLUON PROPAGATOR AND RENORMALIZATION PROCEDURE 



Lattice QCD simulations are performed on a finite 4D torus. Therefore, either the infinite volume limit should 
be taken or the simulations must be performed in a sufficiently large volume. Of course, the precise meaning of a 
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aL (fm) 


3.49 


4.65 






# Conf 


88 


99 







TABLE I: The lattice setup. For the conversion to physical units we took the lattice spacing measure from the string tension 
[8^| . The first set of configurations, i.e. those with /3 = 5.7, were generated by the Berlin-Moscow- Adelaide group and the 
results published in [Q. Note that in their paper, the lattice spacing was taken from tq. The Berlin- Moscow- Adelaide data 
was rescaled appropriately to follow our conventions. 



sufficiently large volume depends on the problem to be addressed. Considering pure Yang-Mills theory and taking 
the mass of the lightest glueball, i.e. M g i ue b a u ~ 1.7 GeV [H, |H| as a typical hadronic scale, the corresponding 
length scale is L ~ 0.1 fm. However, the scale at which nonperturbative physics sets in is already at *~ 1 fm. To 
investigate nonperturbative physics and, in particular, the infrared gluon propagator, one should consider volumes 
well above the 1 fm scale. 

In this work we will analyze the propagator computed from the lattices described in TABLEfl] Of the three (3 values, 
P = 6.0 will be used to perform an extrapolation to the infinite volume limit, whilst the Berlin-Moscow- Adelaide at 
f3 = 5.7 and /3 = 6.2 will be used to cross-check the final results. In what concerns the computation of the gluon 
propagator, we will use standard definitions which can be found in, for example, [l|, iL and as such, it will not be 
repeated here. The gauge configurations were generated using version 6 of the MILC 88j . 

The lattice data for the propagators computed at j3 = 6.0 and j3 — 6.2 were chosen as follows. For momenta higher 
than ~ 1 GeV, only those momenta which survive the conic cut [l[ are used. In this way we avoid the problems 
associated with the breaking of rotational invariance. For momenta below ~ 1 GeV, all momenta were included in 
the analysis. In this way, we hope to have obtained a decent description in the infrared. 

Our simulations are done at different lattice spacings. Therefore, in order to compare the propagators computed 
at different /3 values, the data has to be renormalized. In practice, we have renormalized the gluon propagator, after 
performing a conic cut, by fitting the data to 



D Lat {p 2 ) = Z 



In 



A 2 



P 2 



(11) 



Using the lowest order /3-function and the coefficient of the lowest order anomalous gluon dimension for pure SU(3) 
Yang-Mills theory, extracted from e.g. [8ll l82{. gives 7 = 13/22. The fits to equation (fTTj) were performed for a wide 
interval of momenta \p m i n ,Pmax\- For each lattice, the fitting range was chosen so that x 2 /d.o.f. ~ 1, while keeping 
the largest possible fitting interval. For the various lattices, the fitting range and the quality of the fit can be found 
in TABLE [II] 

The renormalized gluon propagator, 

D(p 2 ) = Z R D Lat (p 2 ), (12) 
is related to the bare lattice propagator Di a t(p 2 ) by requiring that 

D(P 2 )L =II2 = 4, (13) 
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48 


64 








Pmin (GeV) 


2.121 


1.591 








Pmax (GeV) 


5.286 


5.110 








X 2 /d-o.f. 


0.95 


1.02 








Zr 


0.0743(72) 


0.0740(27) 









TABLE II: Ultraviolet fits to equation (H). Note that for /3 = 5.7 and the lattice 72 4 , the x 2 /d.o.f. was never below 3. The 
renormalization constants Zr were computed as described in the text and using \i = 3 GeV as renormalization scale. The 
errors on Zr were computed assuming Gaussian error propagation. 



which defines (part of) a particular momentum subtraction (MOM) scheme. This condition defines the renormal- 
ization constant Zr. Here we chose \i — 3 GeV. The values of Zr are reported in TABLE [TTJ The renormalized 
propagator can be seen in FIG. Q] for j3 = 5.7 and in FIG. [2] for the other j3 values. 



Renormalized Gluon Propagator - 1^ = 3 GeV 




P [GeV] 



FIG. 1: Renormalized gluon propagator for /3 = 5.7 simulations. 



The p = 5.7 data seems to define a unique curve. In this sense, one can claim that finite volume effects are under 
control. On the other hand, the propagators computed with larger j3 values show a small dependence on the volume, 
especially in the infrared region. Note that, for the largest two lattices, despite the larger statistics, the /? = 5.7 data 
displays a kind of "wiggling" structure. It is unclear whether this structure is of any relevance, as the /3 = 6.0 and 
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Renormalized Gluon Propagator - [i = 3 GeV 
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FIG. 2: Renormalized gluon propagator for f3 — 6.0 and /3 — 6.2 simulations. 



P = 6.2 data shows no such fluctuations in the infrared. Of course, fluctuations in D(p 2 ), even if they are small, can 
compromise the quality of the fittings. We also observe that the f3 = 5.7 data is below the ft = 6.0 and j3 = 6.2 data 
for momenta smaller than ~ 400 MeV, as it is clear from FIG. [1] and FIG. 



III. THE LATTICE GLUON PROPAGATOR AND THE REFINED GRIBOV-ZWANZIGER APPROACH 



A. Preliminaries 

In [23|, [24| the gluon propagator, among other things, was investigated by exploiting the Refined Gribov-Zwanziger 
action, and the tree level result (TTU1) was derived. This propagator counts three mass scales: (1) M 2 , related to the 
condensation of the new fields {^ b ,^f\^ b ,<^''} introduced by Zwanziger [13, EH to localize the Gribov-Zwanziger 
action which was nonlocal after the first step of the construction, (2) m 2 related to the (A 2 ) condensate, and (3) 7 4 
multiplying the horizon function, which is introduced to suppress the Gribov copies in the functional integration. The 
differences between the (Refined) Gribov-Zwanziger action and the usual Faddeev-Popov action can only show up in 
the infrared, as the difference between both is of a soft nature, i.e. proportional to j 2 , see the action (jSJ). Indeed, if 
we formally set j 2 = 0, the GZ action reduces to the usual Faddeev-Popov action as the {vjf iV^f, w° h ,o7° fc } fields 
can be integrated out to form a unity. 

In principle, one can expect that the gluon propagator (|10[) should be able to reproduce the lattice data in a certain 
momentum region. Note that being a tree level result, it does not include the logarithmic correction which has been 
observed, for example, when fitting the ultraviolet region. Remember that the logarithmic dependence was explored 
to renormalize the lattice gluon propagator. Therefore, assuming that (fT0|) describes the lattice data, one can expect 
that it will not reproduce the ultraviolet data, the difference being caused by the lack of the perturbative logarithmic 
correction. Anyway, one can explore the lattice results to check if (|10p can reproduce the propagators reported in 
FIGS. [T]and [2] up to a certain maximum momentum, as in the infrared, we expect that the logarithm will "freeze" 
due the presence of infrared mass scales. Furthermore, given the relation between the different mass scales and the 
condensates, by setting either M 2 or m 2 to zero one can check for the corresponding contribution to nonperturbative 
physics. 



B. Gluon propagator and evidence for the d — 2 gluon condensate {A 2 } 

Let us first consider the condensation of the extra fields {<^jf \ <p^ b , uj^ b , wj 6 } . Given that the lattice gluon propagator 
does not vanish at zero momentum, one must have M 2 ^ 0. Indeed, the motivation to introduce the condensate 
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associated with the new ghost-type fields was precisely to be able to have a D(0) ^ [23 
following correspondence between the tree level gluon mass m 2 and the condensate (A 2 ) [2j 



(g 2 A 2 ) 



-Com 2 , 



Co = 



9 N — 1 
13 TV 



We shall need the 



(14) 



which follows from the construction of [65| . From this relation, it is clear that the presence of the condensate requires 
a nonvanishing m 2 . This can be tested setting m 2 = in the tree level expression (|10p and trying to fit the lattice 
data to 



D(p 2 ) 



P 



M 2 



M 2 p 2 + 2g 2 N 1 i 



(15) 



Despite the similar structure of (|10p and (fl"5"j) . the lattice data distinguishes quite clearly the two functional forms. 
Indeed, while (jTOj) is able to reproduce the lattice propagator on a wide range of momentum starting at GeV and 
going up to 1 — 1.5 GeV, in the sense that the corresponding fit have x 2 /d.o.f. < 2, the fits corresponding to (fT5")) 
always have a x 2 /d.o.f . larger than three, and can as such be rejected. 

As an example of a fit with m ^ 0, in FIG. [3] we show the renormalized gluon propagator computed using the 
/? = 6.0 and 64 4 lattice and the fits corresponding to (JTUJ). Although the fits use only the momentum in [0,p ma x], in 
FIG. [3]we show the propagator if one uses (| 101) over the entire momentum region. There is a small difference between 
the lattice data and the prediction of (fit))) in the ultraviolet region which is clearly seen in the gluon dressing function 
- see FIG.SJ As discussed previously, the small observed differences 2 are expected as (fl"0)) does not take into account 
the perturbative logarithmic corrections. 

Renormalized Gluon Propagator - (J, = 3 GeV 




• p= 
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- 64 4 




Fit 




( = 0.948 GeV [x 2 /dof = 1.1] 




- Fit 


" Pma. 


t = 1.243 GeV [Jf/dof = 1.71 





M" = 2.579(60) GeV 
M 2 + m 2 = 0.536(23) GeV 2 
X* = 0.2828(52) GeV"" 

0.948 GeV 



M~ 
M 2 
X. 4 = 



-- 2.448(44) GeV" 
Km 2 = 0.491(18) GeV 2 
0.2708(37) GeV"" 

p = 1.243 GeV 



r l tea t & ±ttb*iUtt£t± 



FIG. 3: Gluon propagator and fit to (|10[) using the momentum range [0,p max ]. Pmax = 1.243 GeV is the largest fitting range 
which has a \ /d.o.f. < 2. The fi gure includes the outcome of the fits for the two fitting ranges considered. 



Our interpretation of the fits to (|T0|) and ([15]) is that the lattice data points towards a nonvanishing gluon condensate 
(A 2 ). In Section III.D, we shall discuss this in more detail and extract an estimate for (A 2 ). 



C. Measuring the scales in the Refined Gribov-Zwanziger gluon propagator using the lattice data 

In this section we aim to investigate the compatibility of the tree level gluon propagator computed using the 
Refined Gribov-Zwanziger action and the lattice data. In particular, we would like to measure the different mass 



2 For the highest lattice momenta p = 7.76 GeV, the measured propagator is 0.01205(32) GeV -2 , while ([TOj predicts 0.0172 GeV -2 . 
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Renormalized Gluon Dressing Function - (j, = 3 GeV 




FIG. 4: The same as in FIG.[3]but for the gluon dressing function p 2 D(p 2 ). The dressing function provides a clear picture of 
the differences between (| 10p and the lattice data in the ultraviolet region. 



parameters in ([TU]) . As discussed previously, it is not expected that (fTU|) is able to describe the lattice propagator 
for the full range of momenta. Therefore, we perform a sliding window analysis, i.e. we shall fit the propagator 
using momenta in [0,p], with increasing values for p. Then, the \ 2 /d.o.f. can be used to establish a maximum range 
of momenta described by (fTO)) - see FIG. [5] For the largest two j3 values and for the largest lattices, the Refined 
Gribov-Zwanziger tree level propagator is able to describe the lattice data well above 1 GeV. In particular, for 
the largest volume, being the (3 = 6.0 and 80 4 case, the lattice gluon propagator can be fitted by (flUf beyond 1.5 
GeV. We draw the reader's attention by noticing that for this particular set of data the "perturbatively" inspired 
expression (TlTj) describes the lattice data starting from 1.5 GeV (see TABLE Hi) . For the smaller /3 = 5.7 simulations, 
the situation is similar, with the exception of the largest two lattices (88 4 and 96 4 ). In the Appendix, we have spent 
a few words about these latter two lattices, and we motivate why we shall keep them out of our analysis. We shall 
however use the other j3 = 5.7 data to check our results later on. 

In FIG[6] we report the result of fitting (|10j) to the renormalized gluon propagator computed from the /3 = 6.0 
and 64 4 lattice data as a function of the fitting range [0,p ma J. Similar plots can be shown for the remaining fits. 
As FIGE] shows, the estimated values for M 2 , M 2 + m 2 and A 4 = 2g 2 Nj 4 + M 2 m 2 are stable against a change on 
Pmax- For each simulation, as a set of values, we choose those which correspond to the largest fitting range with a 
X 2 /d.o.f. ~ 1. For example, for the (3 = 6.0 and 64 4 data, we take p max = 0.929 GeV and M 2 = 2.589 ± 0.068 GeV 2 , 
M 2 + m 2 = 0.539 ± 0.025 GeV 2 , A 4 = 0.2837 ± 0.0059 for a X 2 /d.o.f. = 1.07. When the x 2 / d.o.f. never crosses 
or become to close to 1, such as happens in the smallest fitting lattice volume, we choose the set of values which 
minimizes \ 2 /d.o.f . for the largest possible fitting range. 

In TABLE lllll we report the estimates of the different parameters defining the Refined Gribov-Zwanziger tree level 
gluon propagator for each lattice simulation. The values are plotted in FIG. [7] as a function of the inverse of the 
lattice length L. The data shows a small dependence on 1/L, especially for M 2 + m 2 , and on the lattice spacing, 
i.e. on /3. Nevertheless, for /3 = 6.0, the four volumes can be combined to perform a linear extrapolation to the 
infinite volume limit. 

As a function of 1/L, M 2 is reasonably well described by a linear function. Indeed, the ~\ 2 /d.o.f . of the fit is 2.13, 
giving 

M 2 — 2.15 ± 0.13 GeV 2 , (16) 
in good agreement with the value computed from the largest /3 = 5.7 volume. 
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FIG. 6: Evolution of the fitting parameters with p max for /3 = 6.0 and 64 4 data. 



For M 2 + m 2 , the linear fit gives an infinite volume value of 



M 2 + to 2 = 0.337 ± 0.047 GeV 2 , 



for a x 2 /d.o.f. = 2.04. 



(17) 



For A 4 , the linear extrapolation has a \ 2 /d.o.f . larger than 3. Fortunately, it seems that A 4 shows the smallest 
dependence on 1 / L and the lattice spacing, with the largest volumes providing numbers which are compatible, within 
one standard deviation. Therefore, given the results reported in TABLE IIIII for the largest volumes, one can claim 
that 



A 4 = 0.26 GeV 4 , 



(18) 
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L 


Pmax 


M 2 


M 2 + m 2 


A 4 


xVd.o.f. 








P = 5.7 






64 


1.255 


2.132 ± 0.052 


0.364 ± 0.020 


0.2553 ±0.0051 


0.99 


72 


0.814 


2.017 ±0.097 


0.302 ± 0.028 


0.245 ± 0.011 


1.21 


80 


1.089 


2.151 ±0.047 


0.359 ±0.016 


0.2604 ± 0.0049 


1.55 








P = 6.0 






32 


1.072 


2.82 ± 0.13 


0.652 ± 0.054 


0.2708 ± 0.0096 


1.40 


48 


0.757 


3.07 ± 0.33 


0.71 ±0.10 


0.312 ±0.030 


1.46 


64 


0.929 


2.589 ± 0.068 


0.539 ± 0.025 


0.2837 ± 0.0059 


1.07 


80 


1.103 


2.346 ± 0.043 


0.463 ±0.019 


0.2561 ±0.0030 


1.03 








P = 6.2 






48 


1.419 


2.40 ± 0.11 


0.473 ± 0.045 


0.2677 ± 0.0095 


1.17 


64 


1.476 


2.366 ± 0.066 


0.476 ± 0.027 


0.2721 ± -0.0057 


1.37 



TABLE III: Tree level gluon propagator parameters from fitting the Refined Gribov-Zwanziger propagator (|10[l to the renor- 
malized lattice gluon propagator. The errors reported are statistical and computed assuming Gaussian error propagation. 
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FIG. 7: Parameters for the tree level gluon propagator of the Refined Gribov-Zwanziger action, computed fitting the renor- 
malized gluon propagator, as a function of the inverse of the lattice length L. The reader should remember that, for P = 5.7 
and 72 4 , the lattice data was not well described by the UV fit (|11[) used to define the renormalization constant Zr, see the 
discussion on the renormalization procedure. This can explain the observed fluctuations in the P = 5.7 results. 



which are the reliable digits from the largest two lattices - see the TABLE. The linear extrapolations can be see in 
FIG. [3 We observe that the figures for the j3 = 5.7 data and the linearly extrapolated results are pretty close, giving 
us further confidence in the extrapolation. 

If one uses the extrapolated values, one can write down that 

m 2 = -1.81 ±0.14 GeV 2 . (19) 

Simultaneously we find 



2g 2 N-f i = 4.16 ± 0.38 GeV 4 . 



(20) 
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FIG. 8: The same as FIG.[7]but including the linear extrapolations for M 2 and M 2 +m 2 , which are obtained using the /3 
data. The large volume /3 = 5.7 data serves as a consistency check, as explained before in the text. 



6.0 



Furthermore, assuming that (|10|) describes the infrared gluon propagator, then 

y\/f 2 

D(0) = — = 8.3 ± 0.5 GeV~ 2 . (21) 

The zero momentum gluon propagator computed using the extrapolated values for M 2 and A 4 is in excellent agreement, 
within one standard deviation, with the lattice D(0) computed from lattice QCD for j3 = 5.7 where D(0) ~ 7 — 8.5 
GeV" 2 , P = 6.0 and 80 4 data where D(0) = 8.93 ± 0.47 GeV~ 2 and for f3 = 6.2 and 64 4 data which has a D(0) = 
8.95 ±0.22 GeV" 2 . 

D. Extracting a value for the dimension two gluon condensate {g 2 A 2 } 

In order to obtain an estimate that can be compared with other values available on the market, we shall rely on 
the renormalization group. In particular, we wish to compare with the values ||8} and (|9[) . being 

(g 2 A 2 ) = 5-lltili GeV 2 , (22) 

and 

(g 2 A 2 ) — 4.4 ± 0.4 GeV 2 . (23) 
For the relevant one loop renormalization group equations, we have, in any (massless) renormalization scheme 3 (22TI65T] 

d 2 4 11 N 



2 



V - -w, &>- 3167r 

9 n 2 2 3 TV 

M _ m = 7Q5m , 70 = --^. (24) 



Hence, our estimate ([19]) corresponds to a positive gluon condensate, as using (|14| yields for iV = 3 

{g 2 A 2 } = 3.35 ± 0.26 GeV 2 . (25) 



We recall that the lowest order anomalous dimensions are universal quantities. 
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i.e. a positive gluon condensate. We recall that in this work we have renormalized at a scale [i = 3 GeV. 



The value (|22j) was obtained in the so-called T-scheme, which is kind of MOM scheme compatible with the renor- 
malization prescription (fT"3"|) . at a renormalization scale /i = 10 GeV. The fundamental scale At of this T-scheme is 
related to the conventional MS one through the conversion formula [64] 



Using we have at one loop 



A T = A^e 507 / 792 . 



(26) 



fi—m 



7o 1 



'da"' 2/? ln£ 



-to 



(27) 



Introducing the auxiliary variable £ = In -P- , this can be easily integrated to 



to = m n — = m n 



Co 



In £2- 

At 



(28) 



using the numbers given in (|24|) . The estimate A^jg = 0.224 GeV determined in [64| consequently leads to 



(g 2 A 2 ) 



/j=10GcV 



= 3.03 ± 0.24 GeV 



(29) 



employing (fH|) and (g 2 A 2 ) = 3.29 GeV at /io = 3 GeV as input values. We notice that our estimate is at least in the 
same ballpark as the ones of (|22[) and (|23p . which were obtained in a completely independent way. In these works, it 
was observed that even at relatively large momenta Q 2 , there was a discrepancy between the perturbatively expected 
results, and the lattice estimates for the gluon or ghost propagator and strong running coupling constant. Usually, 
such discrepancies can be accommodated for by nonperturbative power corrections. It was discussed in [62l that 
a power correction proportional to (A 2 ) /Q 2 was necessary to obtain a sensible estimate for e.g. A^g-. In the current 
work, we obtained a lattice estimate for the same condensate (g 2 A 2 ) in a completely different way, hence it is quite 
remarkable a compatible value is retrieved at the end of each analysis. 



E. Extracting an infrared mass scale from the gluon propagator 

As a final effort we would like to estimate an infrared mass scale by using the gluon propagator. A similar attempt 
was done in [8J . The infrared lattice data is well described by equation (fT0|) , which however depends on multiple mass 
scales. For small enough momenta, p ;$ 0.2 GeV, the propagator (fit)]) is well approximated by the so-called pole (or 
Yukawa) propagator fit 

2 _ M 2 _ Z 

^ ' ~ (Af 2 +m 2 )p 2 + A 4 ~ p 1 + m 2 R ' ( ' 



where 



L IR 



M 2 



(31) 



is the infrared mass scale we can associate to infrared pure QCD. Using the infinite volume estimates for A 4 and 
M 2 + m 2 , it follows that 

m m = 771(108) MeV. (32) 

This infrared mass scale is in excellent agreement with the infrared mass scale estimates from large volume SU(3) 
simulations 17] , where a gluon mass in the range 600 - 800 MeV was claimed, and in good agreement with the recent 
value obtained in (86[, where a mm — 651(12) MeV was measured. Furthermore, the value given in equation (|3"Tj) 
agrees also well with the SU(2) result found in (87[, i.e. mm — 856(8) MeV, and it is only slightly larger than the 
SU(2) infrared mass scale derived in Q. In the latter work, it was however noticed that a pole fit like ([3"TJ]) does not 
work out well. Indeed, we observe that in our case, the fitting range is only something like p G [0, 0.2] working in 
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GeV, while the corresponding mass is about 0.85 GeV, so the name of a "pole propagator fit" is a bit misguided of 
course as at p « mjR the fit is already invalid. As an alternative, the authors of [8[ proposed a Gaussian fit in the 
continuum 

D(p 2 ) = Be- {p - po}2/m2 "i , (33) 

identifying from this an infrared mass scale rn? IR . As far as we know, there is no theoretical motivation behind this 
kind of propagator yet. Roughly said, one should identify a mechanism that can generate a momentum-exponential 
into the effective action in the A — A sector. It should also be noticed that in order to write down the expression (|33p . 
an external momentum vector po must be introduced, thereby sacrificing Lorentz invariance . 



IV. CONCLUSIONS 



In this paper, we have shown that 

n( 2] = P 2 + M 2 (M] 

(P ' p 4 + (M 2 + m 2 ) p 2 + 2 5 2 A 7 4 + M 2 m 2 ' 1 ' 

which is the analytical tree level version of the gluon propagator found in the (Refined) Gribov-Zwanziger formalism 



2J |, can describe very well the lattice data for the SU(3) Landau gauge gluon propagator in the infrared. More 



precisely, for momenta up to 1.5 GeV, a good fit was established, which was only possible with nonzero values for 
all mass parameters appearing in (134)) . We discussed their continuum extrapolation, which yielded the following 
estimates: 

M 2 = 2.15 ±0.13 GeV 2 , m 2 = - 1.81 ± 0.14 GeV 2 , 2g 2 Nj 4 = 4.16 ± 0.38 GeV 4 , (35) 
giving the following continuum value for -D(O), 

D(Q) = 8.3 ±0.5 GeV" 2 , (36) 
which is in good agreement with large volume lattice data. 



Since m 2 is related to the (A 2 ) condensate, we were also able to present the value 



{g 2 A 2 f- W G ° V = 3.03 ± 0.24 GeV 2 , (37) 



which compares fairly with other estimates of this d — 2 gluon condensate. 



We conclude that the current work has collected evidence that the Refined Gribov-Zwanziger formalism is perfectly 
well capable of explaining the infrared behaviour of the (lattice) Landau gauge gluon propagator, with its nonvanishing 
zero momentum limit. This is good news, as fitting lattice data is one thing, but one should also be able to explain 
which effects are behind a particular fit. We notice that also certain Schwinger-Dyson results for the same propagator 
describe the lattice data well, see e.g. [U,|35|, perhaps indicative of a close connection between these formalisms and 
the (R)GZ one, a fact already explored in the work [30(. At the same time, we have also provided further evidence 
that one cannot ignore the effects of the dimension two gluon condensate (A 2 ) in the Landau gauge. 



Appendix A: A few words about the /3 = 5.7 data at volumes 88 4 and 96 4 



The gluon propagator computed at j3 = 5.7 and volumes 88 and 96 turns out to be rather problematic to fit. A 
closer look shows that for these largest two lattices the data fluctuates quite strong in the infrared, see FIG. [9] If one 
ignores these data points, the data does behave similarly as for the (3 = 6.0 and (3 — 6.2 simulations. For example, 
for the largest volume, removing the smallest five momenta, i.e. taking into account only p > 214 MeV, the largest 



4 Or more precisely, rotational invariance as we are working in Euclidean space. The breaking is evident as the resulting propagator is no 
longer a function of the invariant p 2 . 
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fitting range with a \ 2 1 d.o.f < 1.6 corresponds to a p Tl 



1.587 GeV. 



The observed infrared fluctuations can be understood from the combination of the (poor) statistics and the 
particular choice of j3. 



In FIG. [TO] we show the results of fitting the 88 4 and 96 4 propagators to equation (fTUj) in the range [0, p ma x] as 
a function of p ma x- The plots only show the fitting parameters with a x 2 / d.o.f. smaller than three; we recall that, 
typically, one considers \ 2 /d.o.f '. below two. As shown in the figure, the values are not stable against a change of 
Pmax, in sharp contrast with the data shown in FIG. [6] Moreover, looking at FIG. [10] we see that the values grow 
with increasing p max and seem to try to approach the typical numbers reported in TABLE lllll We call the reader's 
attention that in FIG. [5] a similar situation happens for the smallest fitting ranges. Indeed, only for p max larger 
than, let us say, ~ 600 MeV, the fitted parameter values start to become stable. In fact, also for all other lattice 
volumes we did consider in the main text, using smaller values of p ma x would give numerical values smaller than 
those reported in TABLE [TTTl while being unstable against variation around the chosen p ma x. 



The inability to fit the 88 4 and 96 4 data over wider momentum ranges can thus be explained by the infrared 
fluctuations. By performing infrared cuts, one could remove these fluctuations and fit the 88 4 and 96 4 propagators up 
to a Pmax well above 1 GeV. However, given that we want to discuss precisely the infrared region, we do not want to 
perform cuts at low momenta. Besides, if one cuts the infrared data for those lattices, then one should also investigate 
its effect for all other lattices. Because of all that, we choose not to include the 88 4 and 96 4 data in the analysis. 

Renormalized Gluon Propagator - (X = 3 GeV 
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FIG. 9: A zoom of the infrared gluon propagator computed at /3 = 5.7. 
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